Instability of 1-loop superstring cosmology 
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A stability analysis is made in the context of the previously discovered non-singular cosmological 
solution from 1-loop corrected superstring effective action. We found that this solution has an 
instability in graviton mode, which is shown to have a close relation to the avoidance of initial 
singularity via energy condition. We also estimate the condition for the breakdown of the background 
solution due to the overdominance of the graviton. 
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Q\ , Tracing back the history of our presently expanding universe, we are naturally led to the era of high temperature 
t-H ■ and large energy density. When the energy density approaches Planck scale, our general relativistic description of the 

universe is no longer considered as valid since quantum gravitational processes are thought to be significant at such 
' a high energy scale. The leading candidate for the satisfactory theory in the very early universe is the superstring 

theory Jl[], and roughly speaking, two approaches are being made in the attempt to construct the consistent history 

of the early universe. 

Pre-big-bang universe model was proposed to realize an inflation H "before" the initial singularity. Due to 
the scale factor duality Q], the equations of motion in this model have two distinct and disconnected branches of 
the solution, one corresponding to an ordinary Friedmann solution, and the other corresponding to the so-called 
super-inflation solution with increasing Hubble parameter. The super-inflationary stage, which is usually put before 
^f~) ' the Friedmann stage and characterized by the negative-power dependence of the scale factor on time, is regarded 
as a sort of inflation that lessens the difficulties appearing in the standard big bang model. Although this model is 
• advantageous in that the inflation is naturally explained from one of the most promising theories at present, many 
' difficulties, besides the serious graceful-exit problem S, have been pointed out 

Another approach in string cosmology is the search for non-singular cosmological models |^|. Initial singularity 
problem is inherent in the standard big bang cosmology |§, and unless we believe the universe was really generated 
by one sudden blow at infinite temperature, it is natural to think that the classical theory was violated in such an 
extreme condition and that singular situation never happened. Also, getting rid of the infinitely high temperature 
y | suggests possibilities to understand the whole history of our universe without knowing the theory of everything. 
5h i Recently, non-singular cosmological model based on 1-loop corrected superstring effective action fli"|| was proposed 
p2| , p^[ . One of the most interesting features of this model is that its solutions include a super- inflationary stage before 
the ordinary Friedmann-likc universe, which is, just like pre-big-bang case, expected to give a natural explanation for 
the inflation. It should be emphasized that the transition in this model is smooth in contrast with the pre-big-bang 
scenario. In order to concentrate on the behavior of the metric-modulus sub-system which plays the essential role in 
the realization of super-inflationary singularity-free solution, a simpler version of this model was proposed by Rizos 
and Tamvakis |l4| , in which they generalized the form of the coupling and analytically examined the conditions for 
the existence of non-singular solutions. In this letter perturbative analysis is made for this metric-modulus system 
and we show that its non-singular solution is geometrically unstable. 

We start with the 1-loop effective action of the heterotic string with orbifold compactification, which is given by 

EMI: 

S = j d^x^g^R-^D^f-^Daf + ^X^-X^)]^}, (1) 
where R and a are the Ricci scalar curvature, the dilaton, and the modulus field, respectively^]. The Gauss-Bonnet 
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1 We will use in this paper the following convention. The signature of the metric: (—,+,+,+), Riemann tensor: R^avp = 
r M Q /3,y + Ricci tensor: R a/ s — R^a^p. Greek indices denote 4 dimensional space-time while Latin indices run from 1 to 3. 
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curvature is 



r>2 tjilvkX td a T>liv d i r>2 / n\ 

U GB — n K-hvk\ — 4 -ft J^tif + -ft j l/j 

and £(<r) is an even function of the modulus field expressed with Dedekind rj function: 

a 00 

£(<t) = -ln[2eV(ie ,T )] = - In 2 - <r + — - 4 V ln(l - e" 2 " 7 ^). (3) 

3 n=l 

This action provides non-singular cosmological solutions driven by the modulus field for positive A2 fll2fl . These non- 
singular solutions experience a super-inflationary phase, which means the expansion of the universe with increasing 
Hubble parameter. Since the non-singular nature of this model depends essentially on the behavior of the modulus 
field ]T^] , the system is well represented by the one scalar field model : 

S = J d±x^-g{ l -R- l -{Dtpf ±£(<p)I& B }. (4) 

Variation of this action gives equations of motion 

G»u = f'^-,, ~ \v' a V-, a ^ v + ^ [H\„p? afi + (■■ r .,C,, - , (5) 

= Y^ R GB- (6) 

Here, iV denotes derivative with respect to tp, is the Einstein tensor and 

H p l av p :— R^avf3 — RotpS^v + Roa/S^p. (7) 

Assuming the homogeneous and isotropic flat metric 

ds 2 = -dt 2 + a 2 S ij dx i dx j \ (8) 
the equations of motion are rewritten using tp and the Hubble parameter H :— a/ a as 

<p 2 = 6H 2 (l-±HQ (9) 

(2H + 5H 2 )(l-^HO + H 2 (l-±O = (10) 
q \ 

<f + 3H(p+—(H + H 2 )H 2 i v =0. (11) 

Dot (') is the derivative with respect to the physical time t. Thus, along with = aH and Aip = ip, we have four 
first order differential equations for four variables a, H, ip and (p, and there is one constraint ([)[). 

As long as we consider the spatially flat universe, the system has two degrees of freedom for initial conditions since 
the scale factor is unimportant except its relative change. Furthermore, because we are only interested in the region 
ip ~ where the energy condition is violated and the initial singularity is avoided jl4| ] , the initial conditions actually 
have only one degree of freedom. Since it is shown in |lj] that (p and H does not change their signs during their 
evolution, we choose (p and H to be positive in order to describe the expanding universe. The action (^) is symmetric 
under the change of <p's sign so that ip can be set to be positive without a loss of generality. We also assume A to be 
positive, since otherwise the system has no singularity- free solutions 

Homogeneous, isotropic and spatially flat solutions leading to Friedmann-like universe (decelerating expansion) in 
the future are shown in fig.l. Since (p is always positive throughout the evolution of the solutions, larger value of p 
means later in time, i.e. time flows from left to right. We here put A = 1, which fixes the time scale. The equations 
of motion are numerically integrated from the future to the past. In fig.l, we fix the initial (future) value of tp and 
take several different initial values of H . The initial values of tp are determined by H and tp through the constraint 
@. We see that there are 2 classes of solutions: singular and non-singular. The singular solutions (a and b in fig. 1) 
lead to an initial singularity as is the case for the usual big-bang universe. The non-singular solutions (c, d and e) are 
free from the initial singularity, and approach a flat space in the infinite past. It is shown in |l4j that all solutions 
that inhabit in the tp < H > quarter plane regularly continue to tp > H > quarter plane. 

Asymptotic solution of the system can be studied by assuming a power-law behavior for H and tp as 
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H-u^tf, (12) 
ip ~ ip + w 2 ln|t|, (13) 

(14) 

in the asymptotic region t — ► ±oo |p^| . Using equations (||)(|ll|) and the asymptotic form of the function £ 

H~sign(^e^, (15) 
we have following asymptotic solutions for t — > oo and t — > — oo regions respectively: 

Aoo : /3 = -1, wi = ^, u; 2 = J|, (16) 
A-oo : = -2, «i > 0, cj 2 = -5. (17) 

The future solution A M is obtained by balancing the terms which do not contain Gauss-Bonnet contribution, and 
therefore describe the situation where the 1-loop effect is negligible. On the contrary A_oo corresponds to the case 
in which the Gauss-Bonnet contribution is dominant, and requires A > 0, which is acceptable for our non-singular 
solution. These asymptotic solutions are uniquely determined under the conditions mentioned above, that is, the 
solutions are singularity-free and H > 0, ip > 0. To sum it up, the non-singular universe we are considering here 
begins with an asymptotically flat space at the past infinity, and experiences super-inflation until ip crosses zero, and 
then enters into a stage of ordinary Friedmann-like expansion where Gauss-Bonnet effect is negligible. 

For the analysis of perturbation we now look into a typical non-singular solution more closely. To characterize the 
equations of state of the system we introduce a new variables T defined by 

p eS =(T-l)p eS . (18) 

Effective energy density p e g and effective pressure p e Q used here are defined by the components of Einstein tensor as 
p e g = —G°o and p e g = \G l i, respectively. Using the background variables we can express T in a more explicit form: 

r = -5§- < 19 > 

Evolutions of H, p and T in a typical non-singular solution are shown in the fig. 2. The friction 3H + a/a is for later 
discussion. In the far future T approaches 2, which indicates that the equation of state for this system in later time 
is nearly that of stiff matter, i.e. free scalar field. On the other hand, F takes large negative values in the far past. 
This peculiar behavior of T has a close relation to the non-singular nature of the solution. By means of this T, weak 
and strong energy conditions p e g- + p e g- > 0, p e g- + 3p e g- > can be written asT>0, T>2/3, respectively. Thus, 
as can be seen from the fig. 2, violation of those energy conditions, which is necessary to avoid the initial singularity, 
occurs near and before the peak of the Hubble parameter H at (p ~ 0. 

We now consider the perturbation of our model (^) to analyze the stability of the non-singular solutions. We include 
only the tensor perturbation here0, and write the full metric as 

ds 2 = -dt 2 + a 2 {5 l] + h^dx'dy 3 . (20) 

For hij we assume the transverse-traceless conditions h l i = 0, hij^ = 0. Equations of motion for hij are obtained 
by perturbing the Einstein- like equation (|^), and non- trivial one comes from the spatial part. Expanding hij with 
transverse-traceless basis tensors as 

hij = h + (k,t)e +i j(k) + h x (k,t)e xij (k), (21) 
we obtain the same equation for each polarization mode 



2 We have also analyzed scalar and vector part of the perturbation using the gauge-invariant method but no instability 
was found. 
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h+(3H + -)h+-^-(l-^i)h = 0, (22) 
a a^a 1 

where a is a function of the background variables defined by 

Equation ( p2| ) describes the evolution of the geometrical fluctuation in our non-singular model. Of course, this is 
nothing but the equation for the gravitational wave in FRW background if A = 0. The background equation (|10| ) and 
the definition of T @ lead to 1 - f £ = (3r - 5) a so that the perturbation equation (E2J) can be expressed as 

a ■ k 2 

h+(3H+-)h+—(3T-5)h = 0, (24) 
a er 

which is suitable to discuss the stability of the system in relation to the energy conditions. In this expression the 
third term determines the stability of the system, and the second term is regarded as a friction. One can see that the 
system is unstable for tensor perturbation when 

r < \ (25) 

and that strong instability appears in small scale. We can combine this result with the discussion about the energy 
conditions to reach the statement, that is, Cosmological solution of this model is geometrically unstable as long as 
it is non-singular. In order to avoid the initial singularity, the strong energy condition T > 2/3 must be violated. 



Then, since F < 2/3 < 5/3, condition (25) is satisfied and the system becomes unstable for tensor perturbation. 
Therefore, the non-singular nature of the background space-time inevitably requires its geometrical instability, which 
is characterized by the condition (|2^). 

Numerical solutions of perturbation for different wave numbers in the flat background are shown in the fig. 3. The 
ordinate is in logarithmic scale, and we can easily recognize that the instability arises in the super-inflationary phase. 
In these calculations the initial value of h is set to unity so that logioh starts from zero. Our statement can be 
confirmed by comparing fig. 2 and fig. 3. One interesting point we can notice is that the effect of the friction term 
3H + a/a is not necessarily negligible. The friction term is negative before the peak of the Hubble parameter and 
this negative friction enhances the instability. After the peak, the friction turns to be positive and it weakens the 
growth of the perturbation. This is why the instability seems to cease before T gets greater than 5/3. 

One might wonder whether this instability continues from the infinite past or not. To answer this question we can 
use the asymptotic background solution which tells us the asymptotic form of T in the infinite past. Substituting 
( |l7| ) into (|l9|), we have the asymptotic form of T as T ~ 4t/3uJi. Since lo\ is positive, T decreases indefinitely in the past 
and the energy conditions can never be recovered. From the condition (|25|), therefore, the geometrical instability never 
disappears even if we go to the infinite past. To obtain the asymptotic solution of h in the past asymptotic region, we 
estimate each term of the equation(p4[). As a/a behaves like l/\t\ and H ~ ui/|t| 2 in the far past, the friction term in 
( pi| ) can be neglected. In the third term, the scale factor a behaves like a = exp J Halt ~ exp(wo — ^r) ~ expcjQ =: io 



(u>0 and cio are constants). Therefore, ( |24j ) becomes 

4fc 2 

h + zth = 0. (26) 

and the asymptotic solution of h can be expressed using Airy functions as h(k,t) ~ C\Ai(t) + C2-Bi(r), where 

V ^1"0 1 1 

The fact that the perturbative instability continues from the infinite past does not necessarily mean the failure of 
the background solution. To see the criterion for the breakdown of the background solution, we assume a state of the 
perturbation at some finite past time t = ti as h(k,ti) = hi(k) and h (k,ti) = h'^k) (prime denotes differentiation 
with respect to r). Then, C\ and C2 are determined and we can write 

h(k, t) = n (Bi'{Ti)hi{k) - Bi( n )ti t (k)) Ai(r) - vr (Ai'( Ti )hi(k) - Ai( Tl )ti t (k)) Bi(r), (27) 



where r, = ?/ 4fc2 2 |^|- In order to keep the perturbation smaller than the back ground metric through the evolution 



LOiao 



of the solution, h(k, t) has to be smaller than unity when t ~ 0, which is sufficient since the growth of the perturbation 
terminates around the Hubble peak. Then, discarding the decaying mode and picking up the dominant term {hi{k) 
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and/i£(fe) are assumed to be comparable in magnitude), the background solution breaks down if and only if h(k, 0) = 
irBi 1 (Ti)hi(k)Ai(0) 3> 0(1). We can use the asymptotic form of the Airy function to put this breaking down condition 
expressed as a condition for the value of the perturbation at time tj 



hi(k) » Tr-^i(O)- 1 ( 12 | ti |-lexpf--^|t i |i) ! (28) 



where Ai(Q)^ 1 w 2.817 ~ 0(1). Therefore, premising the existence of some natural cut-off for k, we can say that 
the background system breaks down due to the overdominant graviton if the amplitude hi(k) of at least one mode k 
satisfies the condition ( |2"§| ) at time ti. 

We have shown that the non-singular cosmological model proposed by Antoniadis, Rizos and Tamvakis is 
unstable for tensor perturbation, and that the energy condition breaking required to obtain the non-singular solutions 
is a sufficient condition for this instability. We also estimated the condition for the break down of the background 
solution due to the overdominant graviton(p8|). 

Let us conclude this letter by giving some comments on our result. First, we would like to comment on the nature 
of the instability found here. One can say that this instability has an opposite nature to Jeans instability, in the sense 
that it appears in small scales rather than in large scales. Apart from its characteristic spectrum, this instability is 
distinctive in that it appears only in the tensor part. Since it involves no growth of scalar field perturbation [Tq| , this 
instability is purely a "geometrical" one. If we assume that the effective action ([!]) is valid in the early stages of our 
universe, and if none of the hi{k) satisfies the condition (^8|) at time ti so that the background solution is unaffected by 
the perturbation, this small scale instability may indicate generation of many small primordial black holes, typically 
of the Planckian scale. Next, although we have treated only the spatially flat case so far, our claim also holds for 
the metric with non-zero spatial curvature ^. However, if the background universe has positive or negative spatial 
curvature, the asymptotic behavior of the system in the infinite past is not the same as in the flat case ]13[ . Therefore, 
if the energy conditions should be recovered in the finite past, the instability might no longer be eternal. Lastly, we 
have to mention the possibility that the effective action we use here (|l|) might not be suitable for the model of the 
early universe. Since the 1-loop correction changes the tree-level solutions drastically, it is sensible to consider that 
the higher loop effect might change the nature of the cosmological solutions from that of 1-loop action. Our results 
of rather pathological feature of the instability might be a consequence of our neglecting higher corrections. To take 
all these loop effects into account, we have to wait for the accomplishment of a non-perturbative string theory. 
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Figure captions 



Fig.l 

Flat(/C = 0) background solutions that lead to decelerating expansion in the future. Numerical integration is per- 
formed from future to past, and the figure shows the phase diagram of Hubble parameter H versus modulus field 
ip. Time flows from left to right since ip is always positive. There are two classes of solutions: singular solutions (a 
and b), and non-singular solutions (c,d,e). One can see that the avoidance of the initial singularity takes place when 
ip ~ 0. For detailed discussions see references [||| and Q. 

Fig.2 

Evolutions of (i)H, (ii)ip, (iii)3iJ + a/a and (iv) T in a typical super-inflationary flat solution. The origin of time t 
is chosen so that <p crosses zero when t = Q. H has a peak at t ~ 0. tp accelerates from the asymptotically flat space, 
and decelerates at ip ~ 0. The friction term 3H + a/a in the perturbation equation of motion evolves from a small 
negative value, increases in its magnitude, and changes its sign at t = 0, and then decreases in the future. T varies 
from a large negative value in the past to 2 in the future. 

Fig.3 

Evolution of the tensor perturbation in a flat background. Note that the tensor perturbation is plotted with logarithmic 
scale. We set a = 1, h(k) = 1 at the onset of the integration (t = —2), and t = is the time when <p crosses zero. One 
can notice that growing modes appear in the initial super-inflationary phase. The smaller the scale of perturbation 
is, the larger the growth rate becomes. The depression of k ~ 20 solution around t ~ 1.25 indicates the oscillatory 
behavior in the Friedmann-like phase. 
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Fig. 2 (ii) 
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Fig. 2 (iv) 
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